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To elucidate the origin of spin structure-driven anomalous Hall effect (AHE) in pyrochlore com- 
pounds, we construct the t2g-orbital kagome lattice model and analyze the anomalous Hall con- 
ductivity (AHC). We reveal that a conduction electron acquires a Berry phase due to the complex 
d-orbital wavefunction in the presence of spin-orbit interaction. This "orbital Aharonov-Bohm (AB) 
effect" produces the AHC that is drastically changed in the presence of non-coUinear spin structure. 
In both ferromagnetic compound Nd2Mo207 and paramagnetic compound Pr2lr207, the AHC given 
by the orbital AB effect totally dominates the spin chirality mechanism, and succeeds in explaining 
the experimental relation between the spin structure and the AHC. Especially, "finite AHC in the 
absence of magnetization" observed in Pr2lr207 can be explained in terms of the orbital mechanism 
by assuming small magnetic order of Ir Sd-electrons. 

PACS numbers; 72.10.-d, 72.80.Ga, 72.25.Ba 
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I. INTRODUCTION 

Recently, theory of intrinsic anomalous Hall effect 
(AHE) in multiband ferromagnetic metals has been de- 
veloped intensively from the original work by Karplus 
and Luttinger (KL) fl|. The anomalous Hall conduc- 
tivity (AHC) a AH = jx/Ey due to intrinsic AHE shows 
the almost material-specific value that is independent of 
the relaxation time. The intrinsic AHE in heavy fermion 
compounds , Fe Q , and Ru-oxides had been stud- 
ied intensively based on realistic multiband models. Also, 
large spin Hall effect (SHE) observed in Pt and other 
paramagnetic transition metals [7] , which is analog to the 
AHE in ferromagnets, is also reproduced well in terms of 
the intrinsic Hall effect The intrinsic AHE and 

SHE in transition metals originate from the Berry phase 
given by the d-orbital angular momentum induced by the 
spin-orbit interaction (SOI), which we call the "orbital 
Aharonov-Bohm (AB) effect" [ll|. 

In particular, AHE due to nontrivial spin structure 
attracts increasing attention, such as Mn oxides (T2j 
and spin glass systems fl3|. The most famous exam- 
ple would be the pyrochlore compound Nd2Mo207 [13- 
\1T\ . Here, Mo 4d electrons are in the ferromagnetic state 
below Tc = 93K, and the tilted ferromagnetic state in 
Fig. [T] is realized by the non-coplanar Nd 4/ magnetic 
order below Tn « 30K, due to the d-f exchange in- 
teraction. Below Tn, the AHC is drastically changed 
by the small change in the tilting angle 9 of Mo spin; 



9^ < 10 in the neutron-diffraction study 16:]. This be- 
havior strongly deviates from the KL-type conventional 
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FIG. 1: Tilted ferromagnetic state in the kagome lattice. Blue 
circles are Mo ions. Arrows at Mo sites are the tilted ferro- 
magnetic exchange field. 



behavior a ah oc oc 1 — ^9^. Moreover, the AHC 
given by the spin chirality mechanism [Til , [20| . which 
is proportional to the solid angle sa • (sb x sq) cx 9^ 
subtended by three spins, is also too small to explain 
experiments. Moreover, (9'^) takes the minimum value 
under H ^ i Tesla according to the neutron-diffraction 
study [iBl, whereas the AHC monotonically decreases 
with H. Thus, the origin of the unconventional AHE 
in Nd2Mo207 had been an open problem for a long time. 

Very recently, this problem was revisited by the present 
authors by consideriiig the d-orbital degree of freedom 
and the atomic SOI [2l|, and found that a drastic spin 
structure-driven AHE emerges due to the orbital AB ef- 
fect, in the presence of non-collinear spin order. Since 
the obtained AHC is linear in 0, it is much larger than 
the spin chirality term for \9\ ^ 1. In Ref. [2l[, we con- 
structed the t2g orbital kagome lattice model based on 
the spinel structure (ArMo204): Although Mo atoms in 
XM02O4 and A"2Mo207 are equivalent in position and 
forms the pyrochlore lattice, positions of O atoms in 
Ar2Mo207 are much complicated. 

In this paper, we construct the t2g kagome lattice tight- 
biding model based on the pyrochlore structure, by tak- 
ing the crystalline electric field into account. We find 
that the orbital AB effect causes large 0-linear AHC, re- 
sulting from the combination of the non-collinear spin 
order (including orders with zero scalar chirality) and 
atomic SOI. The realized AHC is much larger than the 
spin chirality term due to non-coplanar spin order, and 
it explains the salient features of spin structure-driven 
AHE in Nd2Mo207. We also study another pyrochlore 
compound Pr2lr207, and find that the orbital AB effect 
also gives the dominant contribution: We show that im- 
portant features of the miconventional AHE in Pr2lr207, 
such as highly non-monotonic field dependence and resid- 
ual AHC in the absence of magnetization, are well repro- 
duced by the orbital AB effect. 
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TABLE I: Coordinates of Mo and O in pyrochlore structure 
as shown in Fig. Oin the xj/z-coordinate. 



FIG. 2: Pyrochlore structure. Blue (white) circles are Mo (O) 
ions. The Mo ions on the [111] plane form the kagome lattice. 




FIG. 3: Configurations of Moa-d and Oi-g. The xsyszs- 
coordinate is defined by the surrounding O tetrahedron. 



The paper is organized as follows: In Sec. II, we intro- 
duce the pyrochlore-type t2g orbital tight-binding model 
and the Hamiltonian. We give the general expressions 
for the intrinsic AHC in Sec. Ill, and explain the orbital 
Aharonov-Bohm effect in Sec. IV. The numerical results 
for Nd2Mo207 and Pr2lr207 are presented in Sec. V 
and VI, respectively. In Sec. VII, we make comparison 
between theory and experiment. 



II. MODEL AND HAMILTONIAN 

First, we introduce the crystal structure of the py- 
rochlore oxide A2B2O7: It has the face centered cubic 
structure, in which two individual 3-dimensional net- 
works of the corner-sharing A4 and B4 tetrahedron are 
formed. In this paper, we mainly discuss the AHE in 
Nd2Mo207, and Pr2lr207 is also discussed in section IVTl 
Figure [2] represents the Mo ions (Blue circles) and O ions 



Ion 


Site 


Coordinate 


Mo 


A 


(1/4 ,0, 0) 




B 


(0, 1/4, 0) 




c 


(0 1/4) 




D 


(1/4, 1/4, 1/4) 





1 


(1/8, 1/8, -1/16) 




2 


(1/8, -1/16, 1/8) 




3 


(-1/16, 1/8, 1/8) 




4 


(5/16, 1/8, 1/8) 




5 


(1/8, 5/16, 1/8) 




6 


(1/8, 1/8, 5/16) 



(White circles) in the pyrochlore structure. The [111] Mo 
layer forms the kagome lattice. The Mo 4(i-electrons give 
itinerant carriers while the Nd 4/-electrons form local 
moments. 

We construct pyrochlore type t2g-orbital tight binding 
model in the kagome lattice for Mo Ad electrons, where 
the unit cell contains three sites A, B and C in Fig. [Jl 
The coordinates of Mo and O are shown in Table [J 2^] , 
and the quantization axis for the Mo c?-orbital is fixed by 
the surrounding Og octahedron. To describe the c?-orbital 
state, we introduce the (a;?/z)H-coordinate for S = A, B, C 
sites shown by FigO The (a:yz)5;-coordinate is defined 
by the surrounding Og ions. In the case of the {xyz)^- 
coordinate, we choose xa, Va and za axes as Moa — >'0i, 
MoA -^02 and Moa — >'04 direction, respectively, in Fig. 
[31 We also choose the {xyz)B- and (a;yz)c-coordinates in 
the same way. 

Moreover, we introduce the ATFZ-coordinate on the 
kagome layer shown in Fig. |D^b). We choose X axis 
as Mo A — >Mob direction and Y axis is perpendicular 
to X axis on the kagome layer. A vector {nx,ny,nz)E in 
the (a;j/z)2-coordinate is transformed into [nx, ny, nz] in 
the XyZ-coordinate as {nx,'riy,nz)s = [nx,nY ,nz\OE., 
where the coordinate transformation matrix O- is given 

by 



Oa = 



Or = 



On. - 



3V6 



3^/6 



3\/6 



Ays ^/i 

-2 7 1 , (la) 
V-%/2 -\/2 5V^/ 

'-VS -4^/3 -\/3\ 

7 -2 1 , (lb) 

-3V3 SVS \ 

-5 -5 -2 . (Ic) 

-V2 -V2 5\/2/ 



Arrows in Fig. IH^a) represents the local effective mag- 
netic field at Mo sites, which is composed of the ferro- 
magnetic exchange field for Mo 4c?-electrons and the ex- 
change field from Nd 4/ electrons. Under the magnetic 
field parallel to [111] direction below T/y, the direction 
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FIG. 4: (a) Umbrella like locale exchange field at Mo sites 
represented by arrows. A unit cell contains sites A, B, C. (b) 
Kagome lattice, aij = A,B,C) is a half Bravais vector, 
(c) First Brillouin zone in a hexagon shape. 



where c' is a creation operator for 4(i-electron on Mo 
ions while the field h arises from the ordered Nd mo- 
ments, which are treated as a static, classical background. 

{a, (3) and (ct, cr') represent the sites, t2g-orbitals 
and spins, respectively. Hereafter, we denote the t2g- 
orbitals {xy,yz, zx) as (1,2,3) for simplicity. The first 
term in eq. ([2|) describes electrons hopping, tiaj/3 is the 
hopping integrals between (i, a) and (j, /?). The direct d- 
d hopping integrals are given by the Slater-Koster (SK) 
parameters (dda), (ddn) and (ddS) In the present 

model, however, SK parameter table given in Ref. [l^l 
is not available since the d-orbitals at each site are de- 
scribed in the different coordinate as shown in Fig. [S] 
In Appendix A, we will derive the hopping integral be- 
tween the sites with different coordinates. The second 
term in eq. ^ represents the Zeeman term, where hi is 
the local exchange field at site i. fi^ = -~2s is the mag- 
netic moment of an electron. Here, we put fiB=i- The 
third term represents the SOI, where A is the spin-orbit 
coupling constant, and I and s are the d-orbital and spin 
operators, respectively. 

The Hamiltonian in Eq. ([2|) is rewritten in the mo- 
mentum space as 

H = Y,ClHkCk, (3) 
k 

where k summation is over the first Brillouin zone in Fig. 
IH^c), and 

= (^Mt' "^Lii' ^fc,it' ^Lii' ''Lii' ''Lit 

"fe,2t' "fe,2i' "fc,2t' "fc,2i' ^k,2i^ ^k,2t 
'^fe.at' '^fc,34.' ^fc,3t' ^Lsi' ''1,34.' ''Lst-'' ("^^ 

Here and hereafter, we denote the creation operators at 
sites A, B and C as a'l^^cn^k,acr ^-^d c^^^, respectively. 
Hk is given by 18x18 matrix: 



of the local exchange fields at sites A, B and C in the 
XyZ-coordinate are (0a = 7r/6,6'), {(pB = 57r/6,6') and 
{4'C = — 7r/2,6'), respectively. In Nd2Mo207, the tilting 
angle 6 changes from negative to positive as H increases 
from -f Tesla, corresponding to the change in the spin- 
ice state at Nd sites [IGj, i23j . 

Now, we explain the Hamiltonian. The Hamiltonian 
for the t2g-orbital kagome lattice model is given by 

^ = ^io-. jf)c\^ ,jCjli,a — hi ■ [fJ'da.a' cl^^^Cia.a' 

+ ^ 51 t']".^ ■ NT,^'cL,aCj/3,„ (2) 



Hk21 Hk22 Hk23 I J (5) 
HkSl Hk32 Hk33/ 

where Hkap is a 6x6 matrix with respect to (S,(t). 
Here, we divide the Hamiltonian ([3]) into four parts: 

Hkai3 = Hiap + + H^n + H^p^, (6) 

where we added the crystalline electric field potential 
term H^f^ to Eq. ([2]). The kinetic term H^^p is given 

by 
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PAq,C/3 \ 



/3At /3A4, ^Bt ^Bi 

aAt / PAa,B/3 

aA| PAq,B/3 PAafiH 

aB t Pbq,a,9 Pbq,c/J 

aB I PBa.Afi PBa,Cf} 

aC t PCa.Afi PCa,B/3 

aC I V PCa.Afi PCa3/3 

I 



(7) 



where Pia,jp = '^tiajp cos{k ■ aij ) and a.y is a half Bravais 
vector in Fig. HJb). 

The Zeeman term H^g is given by j20j 



/yzc 



Rf).s — 





/A 


0B 


aA 







hoSap aB 


R( 


*,57r/6 


aC 
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; 




t / cose 


sin eC'f' 




I I sine'e'' 


^ -cose 





/3C 



Rb-tt/2 



,(8) 



(9) 



where /io=|/ii|, and is a Kronecker's deha. 

Now, we consider the SOI term iJ'^ in Eq. (jH]). For con- 
venience in calculating the AHC, we take the Z-axis for 
the spin quantization axis. Then, 2s = 2[sxiSytSz] is 
given by the Pauli matrix vector in the XYZ-coordinate. 
>^ 



operator in the (xyz)H-coordinate, which is given by the 
relationship {s^,Sy,s^) — [s x , sy , s z]Oe. and Eqs. pa|) - 
Px|). In the (a;?/z)H-coordinate, the nonzero matrix ele- 
ments of I are given as (3|ia;|l) = (lllu|2) = (2|?z|3) — i 
and their Hermite conjugates [13, [25|. Thus, the ma- 
trix elements {H^ j3)sa,E'a' = (Sacr|if'^|S/3cr') • fe.s' for 
(a, 13) — (3, 1) are given as 

iX 



{H^,l)Acr,Aa' 



3^/6 
i\ 

3\/6 
i\ 



a|(4y3sx-2sy- V2s2)|a'), 

a|(-V3sx + 7sy- V2sz)|a'), 
cr|(-3-y3sx — 5sy — -\/2s2)|(t'). 



To derive the i?^n, however, we have to express the spin 



Thus, the 3-1 component of the third term in Eq. 
becomes 



-^31 



6V6 





lA t 


lAi 


IBt 


ib; 


ict 


ICi 


3A t 


/ -\/2 
















3a; 


4V3- 2i 


V2 














3B t 








-V2 










3Bi 








-%/3-f 7i 


V2 








3Ct 














~V2 


-3\/3 + 5i 


3c; 


I 











-3\/3- 5i 


V2 



The 1-2 and 2-3 components are calculated in a similar way. The obtained results are given by 



Hl2 



iX 
6\/6 



Hi 





2A t 
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2Bt 


2Bi 


2Ct 




2Ci 


lA t 


/ -V2 


















ia; 


V3 + 7i 


V2 
















IBt 








-^/2 


-4Vi + 2i 










ib; 








-4v^- 


2i V2 










ict 














-V2 




3V3 + 


ic; 


V 











3V3- 5i 






3At 


3Ai 


3Bt 


3B ; 3C t 


3c; 






2A t 


/ 5^2 










\ 




2a; 


VS + i 


-5%/2 















2Bt 








5^2 


-V^-i 









2b; 








-V^ + i 


-5^2 









2Ct 











5V2 


2i 






2Ci 
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-2i 


-5^2 / 







5i 



(11) 



(12) 



(13) 
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Finally, we consider the crystalline electric field Hamil- 
tonian H^^'^ , which describes the splitting of t2g level 
into two levels aig (non-degeneracy) and e'g (two-fold de- 
generacy) by the trigonal deformation of MoOg octahe- 
dron. The crystalline electric field Hamiltonian in this 
case is given by 



(14) 



The eigenvalues of TJ^ef g^^j-^ gj^g g^j.g 2E0 for aig 
state; \aig) = -^ilxy) + \yz) + \zx)), and -Eq for e'g 
states; le^) = 7^ (lyz) - \zx)) and je^) = 7g(2|i2/) - 
\yz) — \zx)). Thus, the crystalline electric field splitting 
between aig and e'g is 3|£'o|. 



III. ANOMALOUS HALL CONDUCTIVITY 

In this section, we propose the general expressions for 
the intrinsic AHC based on the linear-response theory. 
The Green function is given by a 18 x 18 matrix: Gk{e) — 
((e -I- /i)l — Hk)~^, where /i is the chemical potential. 
According to the linear response theory, the AHC is given 
by CTAH = ctah + ^ 



^AH A 



C^AH 



TI 
f^AH 



k 

-T 



Attn 



deTv 

fc 

-JxG^'jY^-iR^A) 
ae 



e=0 

dG^ 



(15) 



de 



JyG' 



(16) 



Here, G^'"^\e) = Gfe(e + (— )*7) is the retarded (ad- 
vanced) Green function, where 7(> 0) is the quasiparticle 
damping rate, jkfi = -edHk/dkf^ = -ev^ {fi = X,Y) 
is the charge current, where — e is the electron charge. 
Since all the matrix j^^ is odd with respect to fc, the 
current vertex correction due to local impurities vanishes 
identically [1, Thus, we can safely neglect the cur- 
rent vertex correction in calculating AHC in the present 
model. In the band-diagonal representation, eqs. 
and pop are transformed into 



_ ^ \^ Ami, In 

^^"^ ~ 27rN ^ 



1 



II 
C'AH 



2TrN 



Im 



1 



(e-i?i + i7)2(e-£;™ + z7) 
1 



(18) 



TABLE II: Phases for t2„ orbitals 











H=A,B,C 


$0 1 J. , fig 

^ 26 Y 2 





$0 4- ±, fie 

za: ^ 13 Y 2 



at zero temperature. Here, I and m are the band indices, 
and we dropped the diagonal terms I = m since their con- 
tribution vanishes identically. We perform the numerical 
calculation for the AHC using Eqs. (IT7)) and (|18l) in later 
section. 

cr^jj and CT^jj are called the Fermi surface term and the 
Fermi sea term, respectively. According to Refs. [Tol.[26|. 
cr^jj can be uniquely divided into cr^y and the Berry 
curvature term CTah- The intrinsic AHC is given by cr;^^ 
when 7 — since ct^h + ^ah = 0- general cases, 
however, the total AHC is not simply given by cr^J^ since 
^AH + ^AH finite when 7 7^ or ji/j,n 7^ 1- Therefore, 



we calculate the total AHC ctah = o-ah 
paper. 



'AH 



this 



IV. ORBITAL AHARONOV-BOHM EFFECT 

Before proceeding to the numerical calculation for the 
AHE, we present an intuitive explanation for the un- 
conventional AHE induced by the non-coUinear local ex- 
change field hi. For this purpose, we assume the strong 
coupling limit where the Zeeman energy is much larger 
than the kinetic energy and the SOI [2l[. The t2g energy 
levels are split into the two triply-degenerate states by 
the Zeeman effect, as shown in Fig. [5l Its eigenstate for 
—ho/ 2 are given by 



|a) = sin -|a t) + e""^ cos -|a i). 



(19) 



where a = xy,yz,zx. In addition, we assume the SOI 
is much larger than the kinetic energy. Since fi^ = —2s, 
the SOI term at site i is replaced with (— A/2)Z-ni, where 
Ui = hi/\hi\. Its eigenenergies in the t2g space are and 
±A/2, as shown in Fig. [5] The corresponding eigenstates 
are given by [2l| 



l"±) 



n^lxy) + Hxlyz) + ny\zx), 



(20a) 



/2{nl + 



\vz) 



-{nxUz ± iny)\xy) 



{n^riy =F in^)\zx)], (20b) 



where n ~ {rix, ny, nz)B. in the (a;yz)H-coordinate is given 
by [sin 9 cos 0, sin 9 cos 0, cos 6']Oh. In the complex wave- 
function |n_), the phase of each d-orbital within the 9- 
linear term is given in Table II. 

Here, we explain that the 6'-dependence of the d- 
orbital wavefunction |n ) gives rise to a prominent spin 
structure-driven AHE [21| . Figure [6] shows the motion 
of an electron: (a) moving from |B;yz) to \C;zx), (b) 
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transferring form \C; zx) to \G;yz) at the same site, and 
(c) moving from \C;yz) to \A;zx). Here, we assume 
that the electron is in the eigenstate |ti_) at each site. 
The total orbital phase factor for the triangle path along 
A— 7>B— >-C— >-Ais given by the phase of the following 
amplitude: 



To,b = (A;n_|iJ*|C;n_)(C;n_|iJ*|B;n_ 
x(B;n_|7?*|A;n_), 



(21) 



where is the kinetic term in the Hamiltonian. 
For simplicity, we take only the following largest hop- 
ping t = {A;zx\H*\B;yz) = (B; zx\H^\G;yz) = 
(C; 2x|iJ*|A; j/z), and assume that t is real. Consider- 
ing that |S;n_)(S;n_| 3 |S; za;|i exp(-i($°^ 

T^w|6')) for S=A, B, and C, the hopping amplitude is 



13 Y 2 

expressed as 



orb 



\Torh\e 



(22) 



where $o = 27rfi/|e| is the flux quantum, and 
27r$orb/$o = 3$°^ + (9/13)^3726' is "the effective AB 
phase" induced by the complex d-orbital wavefunction. 
The large 6'-linear term in $orb gives rise to the large 
spin structure-driven AHE in Nd2Mo207. 

Note that (B; n_ | A; n'_ ) is not actually a real num- 
ber if n ^ n', since the rotation of the spin axis induces 
the phase factor; see Eq. (IT^ . This fact gives rise to the 
effective flux due to the spin rotation 4>spin; — 47r$spin is 
given by the solid angle subtended by wa, ^^b and nc 
[igj . Thus, the total flux is given by $tot = $orb + ^^pUy 
spin oc 9"^ is negligible for \9\ <^ 1 



However, $ 

Since all the upward and downward ABC triangles in the 
kagome lattice are penetrated by <&tot ~ ^orb, the orbital 
AB effect induces prominent spin structure-driven AHE 
in Nd2Mo2 07. 



V. NUMERICAL STUDY 

In this section, we perform numerical calculation for 
the AHC using Eqs. (ITT]) and p^ . using realistic model 
parameters. We use two SK parameters between the 
nearest neighbor Mo sites as SK(— 1.0, 0.6, —0.1) and 
SK(— 1.0, 0.4, —0.1) where we represent the set of SK pa- 
rameters as SK{{dda) , {ddir) , {dd6)) . Hereafter, we put 
the unit of energy Kddcr)! = 1, which corresponds to 
2000K in real compound. The spin-orbit coupling con- 
stant for Mo 4d is A = 0.5 To]. The number of electrons 
per unit cell is = 6 (1/3- filling) for Nd2Mo207 since 
the valence of Mo ion is 4-1- . We choose \hi\ to reproduce 
the magnetization of Mo ion 1.3/iB in Nd2Mo207 Tg']. 

Figure [7] shows the total and partial density of states 
(DOS) for SK(-1.0,0.6,-0.1) at (a) Eo = and (b) 
Eo = —2, with the damping rate 7 — 0.05. For Eq — 
{Eo — —2) we set \hi\ = 0.8 {\hi\ = 3.0) to reproduce the 
magnetization of Mo ion 1.3/j,b [3. The crystalline elec- 
tric field splitting for Eq — —2 in Fig. [7] (b) corresponds 



E 



(a) 



]2g 



■tA/2|"'\^ 



FIG. 5: Eigenenergies for t2g electron under the exchange field 
h and the SOI (-A/2)n ■ I; see Ref. [2l|]. 




FIG. 6: Orbital AB phase given by the complex t2g orbital 
wavefunction at site C. The electron acquires the phase dif- 
ference between zx- and j/z-orbitals (orbital AB phase) via 
the movement (a)— 7>(b)— 7>(c). 



to 3|£'o| ^ leV, consistently with the band calculation 
'27']. In both cases, the states \aig \.) and \e'g t) gives 
large partial DOS near the Fermi level. 

Figure [S] shows the 3rd- 12th bands from the lowest. 
Nine bands near the Fermi level(A^ = 6) are composed 
of \aig I) and \e'g f) as understood from Fig. [Zljb). As 
shown in Fig. |9l the band structure and the Fermi surface 
are hardly changed by varying by 3 degrees. 

Here, we present the numerical results of the AHC 
for two SK parameter sets; SK(— 1.0, 0.6, —0.1) and 
SK(-1.0,0.4,-0.1). We set {Eo,\hi\) = (-2,3.0) or 
(0,0.8); each parameter set reproduces the magnetiza- 
tion of Mo ion 1.3/^B [Hi- We also put the damping rate 
7 ~ 0.001 (clean limit) unless otherwise noted. Here- 
after, the unit of the conductivity is /ha, where h is 
the Plank constant and a is the lattice constant. If we 
assume a = 4A, /ha w 10'^O~^cm~^. In the numerical 
study, we use 512^ fc-meshes. 

Figure [TU] shows the obtained AHC for 
SK(-1.0,0.6,-0.1) at 6* = and ±3°, for (a) a 
wide range of fi [N = 4.8 — 12) and (b) a narrow range of 
/i {N = 5.0 — 6.7). Since the present spin structure-driven 
AHE is linear in 0, a very small causes a prominent 
change in the AHC although the Fermi surfaces are 
hardly changed (see Fig. |9]). The ^-dependence of the 
AHC for other SK parameter SK(-1.0, 0.4, -0.1) is 
shown by Fig. [TT] for N = 4.8 — 6.6. A remarkable 
change of the AHC is also caused by small change in 0. 
Therefore, large 0-linear term in the AHC is obtained by 
using general SK parameters. 

The finite AHC at = is nothing but the conven- 
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FIG. 8: 3rd-12th bands from lowest for SK(-1.0, 0.6, -0.1) 
and 61 = 0. 
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FIG. 9: Fermi surface for SK(-1.0,0.6, -0.1), (a) 6^0, and 
(b) e = 3°. 
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FIG. 10: /.i-dependence of AHC for SK(-1.0, 0.6, -0.1), for 
(a) N = 4.8- 12 and (b) iV = 5.0 - 6.7 . 



tional KL-term. However, obtained 6'-linear AHC devi- 
ates from the conventional KL-term that is proportional 
to the magnetization oc 9^. We stress that the large 
0-linear term in Figs. 1101 and fTT] cannot be simply mrder- 
stood as the movement of Dirac points (or band crossing 
points) across the Fermi level, since the change in the 
band structure hy 9 = ±3° is very tiny as recognized in 
Fig. ini Thus, the origin of the 0-linear term should be 
ascribed to the orbital AB phase [2l| discussed in Sec. 
IV. 



8 



N=5.0 N=6.0 N=6.7 



o 

CO 


/\ 
' \ 
/ \ 


SK(-1.0, 0.4, -0.1). 
Eo=-2 
7=0.001 


^ 
o 
I 
< 








9=3° ■ 


-2 




6=0° - 






— e=-3° 



-1.5 -1.4 -1.3 -1.2 



FIG. 11: ^-dependence of AHC for SK(-1.0, 0.4, -0.1) and 
= 4.8 - 6.6 . 
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FIG. 12: 7-dependence of AHCs for SK(-1.0, 0.4, -0.1) and 
A = 6. The AHC starts to decrease for 7 ~ 0.02, which 
corresponds to p = 0.26mr2cm in the present parameters. 
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FIG. 13: e-dependence of AHCs for SK(-1.0, 0.6, -0.1), A = 
6 and 7 = 0.001, in cases of (-Bo, A) = (-2,0.5), (-Bo,A) = 
(0,0.5), and (So.A) = (-2,0). 
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Next, we discuss the 7-dependence of the AHC. As 
7 increases from 0.001, spike-Uke fine structure in Figs. 
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[TUl and [TT] becomes moderate as recognized in Ref. 
Moreover, the intrinsic AHC shows a crossover behavior, 
that is, the AHC starts to decrease when 7 exceeds the 
band-sphtting A, proved by using tight-binding models 
[2I] or local orbitals approach |28l]. Figure [T2] shows 
the 7-dependence of the AHC in the present model. Line 
(i) represents the total AHC for = 0; o'ah(6' = 0), 
and line (ii) represents the variation of the AHC from 
e = -3° to 3°; |crAH(3°) - crAH(-3°)|. We also calculate 
the AHC for A = 0, which represents the spin chirality 
driven AHC (jf^ ■ Note that o^^^iff) is an even function 
of 61, and uf^ In Fig. [H we plot |2ct7j^"(3°)| 
as line (iii). The variation of the AHC from B — —i° 
to 3° due to the orbital mechanism is 100 times larger 
than the spin chirality term in the clean limit. Note that 
the intrinsic AHC follows an approximate scaling relation 
CTAH oc yO^ [1, H nil in the "high-resistivity regime" . In 
Fig. [121 we see that |(tah(3°) — ctahC— 3°)| also follows 
the relation similarly. 



FIG. 14: e-dependence of AHCs for !6l| < 3° at 7 = 0.001. 
The spin chirality term (A = 0) is very small. 



Next, we analyze the overall 6'-dependence of the AHC, 
by ignoring the experimental condition \Q\ <^ 1. Figure 
[T^ shows the AHCs as functions of 9. Solid and dashed 
lines represent the AHCs for Eq = and —2, respectively. 
They have large 6'-linear terms for 6 0, and they take 
finite values even ii 9 = ±7r/2 (coplanar order). Note 
that obtained 6'-dependence of the AHC is insensitive to 
the value of Eq. The AHCs for 6* = corresponds to 
the conventional KL-type AHE. Dotted line in Fig. [T3] 
shows the AHC for A = 0, which gives the spin chiral- 
ity term cr^y\ It is proportional to 9^ for small 9, and 
becomes zero when 9 = ±tt/2. Finally, we analyze the 
6'-dependence of the AHC more in detail for \9\ < 3° in 
Fig. m In the case of SK(-1, 0.4, -0.1), the AHC for 
A — 0.5 changes the sign at 6* ~ 2° due to the orbital AB 
effect, and it is more that 100 times larger than the AHC 
for the spin chirality term (A = 0). 
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VI. PralraOr 

In the previous section, we discussed the unconven- 
tional AHE in the pyrochlore Nd2Mo207. Here, we 
discuss other pyrochlore Pr2lr207. Unlike Mo 4d elec- 
trons in Nd2Mo207, Ir 5d electrons are in the paramag- 
netic state. Below 0w = 1-7K, localized Pr 4/ electrons 
form non-coplanar spin-ice magnetic order. Under the 
magnetic field along [111], the non-coplanar structure 
of Pr Ising moments are expected to change from "2in 
2out"(i7 - O.TTesla) to "Sin lout"(i7 > O.TTesla). The 
AHC increases in proportion to the magnetization with 
field from to 0.7 Tesla, whereas it rapidly decreases as 
the spins of Pr tetrahedron change from "2in 2out" to 
"Sin lout" for H > O.TTesla. 

On Ir sites in Pr2lr207, the tilted ferromagnet state 
shown in Fig. [T] is also realized. In Pr2lr207, however, 
the ferromagnetic exchange interaction is absent, and the 
local exchange field hi on Ir ion is composed of only the 
exchange field from the Pr moment; ~ Jdf- Since hi is 
parallel to the sum of the nearest Pr momenta, 9 of Ir 
spin is much larger than the 9 of Mo spin in Nd2Mo207. 
Therefore, the tilted ferromagnetic state with large 9 and 
small \hi\ is realized in Pr2lr2 07. 

Now, we explain the local exchange field on Ir sites 
given by Pr tetrahedron. Details of the derivation of 
these local exchange field are presented in Appendix B. 
In the strong magnetic field along [111] O.TTesla), 
the spins of Pr tetrahedron have "Sin lout" structure, 
and the realized local exchange fields at Ir sites are 
((/)A, (/'B, </'c) = (-57r/6, -7r/6, 7r/2) and 9 = 29.5° in Fig. 
[T5] We denote this Ir spin structure as [3 | 1 t]- In this 
section, we promise that < 9 < tt and — tt < (j) < tt. 
In the intermediate field (~ O.TTesla), the spin of Pr 
tetrahedron can take three types of "2in 2out" struc- 
tures with negative Zeeman energy. If we take one "2in 
2out" structure among three, the exchange fields at Ir 
sites are {(l)A,(f>B,<Pc) = (27r/3, tt/S, 7r/2), 9a,b = 58.5°, 
and 9c = 29.5° in Fig. [151 We denote this Ir spin struc- 
ture as [2 I 2 t]- In real compounds, domain structures 
of three "2in 2out" structures are expected to be formed, 
and the total magnetization is parallel to Z-axis. The 
total AHC will be insensitive to the domain structure 
since ctah's due to three [2 4- 2 f] structures are equiva- 
lent. If we take average of three "2in 2out" structures, 
the local exchange fields belongs to the 120°-structure 
with 9 — 14.4°, as shown in Fig. [T5| We denote this Ir 
spin structure as [2 | 2 f]- As a result, the Ir spin struc- 
ture changes as [2 J, 2 t] (or [2 | 2 f]) ->[3 | 1 t] with 
increasing the field from ^O.T Tesla gradually. 

Here, we perform the numerical calculation for 
Pr2lr207. We put the atomic SOI as A = SOOOK^^ which 
is slightly smaller than the atomic value for Ir [10]. The 
number of electrons per unit cell is = 15 for Pr2lr207. 
We set \hi\ = 20K since \hi\ ~ Jjt is estimated to be 
larger than 14K experimentally [2^. We also put the 
damping rate 7 — 0.001 (clean limit). Figure [THj (a) 
shows the AHC in Pr2lr207 with SK(-1.0, 0.6, -0.4). 



H~o.n 



H»o.n 



Pr 




2in 2out 



Ir 




averaged one 2in2out 



120 -structure 

e = 14.4° 

[2T2T] 




30°-structure 

fl = 58.5°, 29.5° 

[2 i 2 T] 




120"-structure 

B = 29.5° 

[3 i 1 T] 



FIG. 15: The spin configurations in Pr tetrahedron are shown 
in the first line. Three "2in 2out" states are realized under 
H < 0.7 Tesla, and one "Sin lout" state is realized under the 
higher field. These Pr spin configurations induce the local 
exchange fields at Ir sites as shown in the last line. In the 
120°-structure, (0a,(^b,<^c) ~ (— 57r/6, — 7r/6, 7r/2). In the 
30°-structure, {(j>A, (I)b, 4>c) = (27r/3, 7r/3, 7r/2). 



Each line represents the AHC for the 120°-structure in 
Fig. [15] The line with "A = 0" represents the spin 
chirality term. In the case of [2 4, 2 t]i the AHC in 
the present model is 10 times larger than the AHC for 
A = 0. Thus, the orbital AB effect dominates the chiral- 
ity mechanism. The variation of the AHC for [2 | 2 f] (or 
[2 2 t]) — >[S I 1 t] can explain the experimental results, 
ignoring the sign of the AHC. For example, the sign of 
the AHC is changed if Jdf is negative. 

In Fig. [H (b), we put SK(-1.0, 0.6, -x) with x = 
O.S,0.4 and 0.5. Although the KL term at 6* = de- 
creases from negative to positive with x, the overall 9- 
dependence of the AHC is not very sensitive to x. 

Recently, Ref. ^ reports that the AHC in Pr2lr207 
shows a hysteresis behavior under the magnetic field be- 
low w IK. That is, the AHC shows the "residual 
AHE with zero magnetization" in Pr2lr207. In terms 
of the spin chirality mechanism, the authors claimed the 
existence of a long-period magnetic (or chirality) order of 
Pr sites with 12 original unit cells [SO]- However, there is 
no theoretical justification for this complex state. Even if 
it is justified, the origin of the hysteresis behavior is un- 
clear. In addition, the magnetic susceptibility ~ Xpr 
does not show anomaly at experimentally. 

Here, we propose an alternative explanation for the 
residual AHE based on the orbital AB effect: In 
Ln2lr207 with _Ln=Nd, Sm, and Eu, the Ir 5d-electrons 
show magnetic order at =S6 K, IIT K, and 120 K, 
respectively [3l|. Thus, monotonically decreases as 
the radius of Ln ion increases. Since Pr is on the left- 
hand-side of Nd in the periodic table, one may expect 
a finite T^{^ IK) in Pr2lr207. We stress that small 
amount of impurities could induce the magnetic order 
in the vicinity of magnetic quantum-critical-point i32j l. 
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under the constraint '^i ~ ^■ 

Under the exchange field by one of '2in 2out" Pr order, 
the obtained Ir spin structure for J \hi\ is shown in 
Fig. [TH For J > \hi\ under H = Tesla, the Ir spin 
structure is changed to the 120°-structure with 9 = 70.5° 
in Fig. [15] (c), which we denote [2 | 2 t]H->+o structure. 
The obtained AHC under this spin structure is —0.006, 
as denoted in Fig. [15] (a): The experimental residual 
AHC is smaller, since the Ir ordered moment is expected 
to be smaller experimentally. The AHC is reversed under 
H = —0 Tesla since the Ir spin structure is reversed. As a 
result, we can naturally explain the "hysteresis behavior 
of the AHC" below - 1 K reported in Ref. 
We stress that the spin chirality term is zero under the 
[2 I 2 t]ff->-±o structure, since Sa ■ [Sb x Sc) = 0. 

Also, the Ir spins under the exchange field by the av- 
eraged "2in 2out" Pr order for J ^ \hi\ show the 120°- 
structure in Fig. [15] with 6 = 70.5°: We denote this 
structure as [2 | 2 t]_ff^+o- The total magnetization is 
zero since the Ir spin on the apical site (not shown) is 
antiparallel to the Z-axis. In this case, we can also ex- 
plain the "hysteresis behavior of the AHC" below ^ 1 
K. However, the sign of the A HC for [2 | 2 tlff^+o is dif- 
ferent from that for [2 | 2 under the positive H. 



VII. DISCUSSION 



A. Comparison between theory with experiments 



6»A = 90° (?ri = 90° 

[2 I 2 T]h->+o 

FIG. 16: e-dependence of AHC in PralraOr for (a) 
SK(-1.0, 0.6, -0.4) and (b) SK(-1.0, 0.6, -x) where x = 
0.3 ~ 0.5. Each line represents the AHC for the 120°- 
structure in Fig. 1151 (c) Ir spin structure below under the 
weak exchange field hi due to "2in 2out" Pr spin-ice order; 
[2 J, 2 tjir-f+o- The total magnetization of the Ir tetrahedron 
is zero since Sc = —Sb \\ CD, where site D is the apical Ir 
site. Note that [2 J, 2 ^Jh-^-o spin structure is the reverse of 
[2 i 2 t]«^+o. 

Here, we analyze the Ir spin structure below , consid- 
ering the classical Heisenberg model for Ir tetrahedron 
under the exchange field hi by Pr spins (see in Appendix 
B): 

4 4 

E = J{Y,S^)^-Y.^^■^^^ (23) 

1=1 i=l 

where 5'^ is the i-th Ir spin, and the positive J (~ T^) is 
the antiferromagnetic interaction between Ir spins. When 
J <^ \hi\, then Si is parallel to hi. When J ^ \hi\, we 
have to find the spin configuration to minimize eq. (j23p 



First, we compare the theory with experiments for 
Nd2Mo207 [iME 111 in detail. Under -ff||[lll] be- 
low Tn, ctah monotonically decreases with H from 
OTesla {9 w -1.5°) to 6Tcsla {9 « 1.5°). This mono- 
tonic decreasing in AHC can be explained by the 9- 
linear term in the present model. The relation ~ 
AnRsM^° + ATTR'^Mz'^ describes the experimental results 
weU, where Mf°{M^'^) and Rs{Rs) are the magnetiza- 
tion and anomalous Hall coefficients for the Mo(Nd) mo- 
ment [l3 - [l^ . In this equation, the first term represents 
the conventional AHE that is recognized as the KL mech- 
anism. In contrast, the second term is highly unusual in 
that Nd electrons are totally localized; it represents the 
unconventional AHE due to the non-coUinear spin config- 
uration. As the magnetic field increases from Tesla to 
6 Tesla, M^'^ increases from negative to positive. Since 
9 oc M^'^ [2ll|. the second term corresponds to the 9- 
linear term given by the orbital AB effect. Moreover, the 
AHC for [2 I 2 t]_f/^+o in Fig- [IS] is finite, irrespective 
of the absence of magnetization. 

Next, we compare the present theory with experiments 
for Pr2lr207 [23|. Under the magnetic field along [111], 
the observed AHC increases in proportion to the magne- 
tization with field from Tesla, whereas it decrease with 
H above 0.7 Tesla as the spins of Pr tetrahedron start to 
change from "2in 2out" to "Sin lout" . The peak value of 
AHC around 0.7 Tesla is 17f7~icm-i. The AHCs in Fig. 
\W\ (a) are -llfl-'^cm-^ for [2 ; 2 t], and it is doubled if 
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(a) Orbital mechanism 

ho-S+ 



(b) Spin chirality mechanism 

ho+s- 




ho-S+ 



FIG. 17: Examples of the second-order diagrams for spin 
structure-driven AHCs: (a) orbital mechanism and (b) spin 
chirality mechanism. The conventional AHE (KL-term) is 
given by the diagrams (a) by replacing Xl+S- and fto-s+ with 
XlzSz and /lozSz, respectively. 



Therefore, (t™^ is about 100 times larger than cr^^" in 
Nd2Mo207. This result is recognized in the present nu- 
meraical calculation in Fig. [T3] and in Ref. [21| . 

In Pr2lr207, the relation \d\ ^ 0(1) is realized. Thus, 
we obtain a^'^^/a"^^ ~ ho/X - O{10^) since ho ~ lOK 
and A ^ 3000K. In the present numerical study, how- 
ever, a AH is only 10 ^ 20 times larger than cr^^", as 
shown in Fig. 1161 This discrepancy originates from the 
higher-order correction of large A on the band-splitting 
A: In fact, <TjI^ starts to decrease for A > 600K since A 
increases with A. 

Finally, we comment that (y^^ is not suppressed by 
large crystalline electric field. Since Z±Sip mixes the states 
\aig, a) and \e'g, —a), a™^ will be large if these two states 
occupy large portion of the DOS at the Fermi level. This 
situation is actually realized the presence of crystalline 
electric field, as shown in Fig. [7] (b). For this reason, 
large 0-linear spin structure-driven AHE is realized for 
Eo = -2. 



we put hi — > 2hi. Thus, the variation of the AHC for 
[2 I 2 t] ^[3 I 1 t] in Fig. [H] (a) can explain the exper- 
imental field dependence. The obtained AHC is mainly 
given by the orbital AB effect, and the spin chirality term 
is too small to reproduce experimental values. 



B. Second-order-perturbation theory for spin 
structure-driven AHCs 

Here, we discuss the spin structure-driven AHC based 
on the second-order-perturbation theory with respect to 
A and hg. The present weak-coupling analysis together 
with the strong-coupling analysis in Sec. IV will provide 
us useful complementary understanding. Since their ex- 
pressions in the present model are too complicated, we 
show only some examples of the the second order dia- 
grams for the spin structure-driven AHCs in Fig. 1171 (a) 

cr™]!^ due to the orbital mechanism, and (b) (J^ah '^^'^ ^'^ 
the spin chirality mechanism. In (a) , the spin of conduc- 
tion electron is fiipped by x, y-components of the Zeeman 
term, ho±s^, and the SOI term, Xl±Sz^, and the obtained 
SHC is a^l^ ~ ho±X/A, where A is the band splitting 
near the Fermi level. We stress that this term vanishes 
when C2 rotational symmetry along Z-axis exists [2l|: 
In the present model, the C2 rotational symmetry of the 
simple kagome lattice in Fig. |l](b) is violated by the exis- 
tence of oxygen atoms. In Fig. [I7](b), the spin is flipped 
by /io±s=F twice, and it is given by ctah" ~ ho±/^- Thus, 
tT™H ^^'^ '''ah" proportional to sin 9 and sin^ 9, respec- 
tively. (We note that the conventional AHE (KL-term) 
is given by replacing Xl+S- and Hq^s^ in Fig. [T7] (a) 
with XlzSz and ho^Sz, and cr^g ^ /iqzA/A.) 

In Nd2Mo207, the relation \9\ ^ 0(10"^) is realized. 
Thus, we obtain (j''^^/cr]l^ 9 since A ~ /lo ~ lOOOK. 



C. Summary 

In summary, we studied the AHE in the pyrochlore 
type t2g-orbital model in the presence of non-coUinear 
magnetic configurations and the crystalline electric field. 
Thanks to the SOI, the complex d-orbital wave function is 
modified by the tilting angle 9, and the resultant orbital 
AB phase gives large 6'-linear AHC. This orbital term, 
cr^'j'^, dominates the AHE in Nd2Mo207 since the spin 
chirality term, cr^^", is proportional to 9^ (^ 1). The 
obtained numerical results are qualitatively equal to the 
resuhs in Ref. [2l|. 

In Pr2lr207, cr^''^ also dominates cr^j^" since the SOI 
for Ir 5(i-electron (A ^ 3000K) is much larger than the 
d-f exchange interaction {Jdf ^ 20K). In particular, the 
present orbital mechanism can explain the "hysteresis be- 
havior of the AHC" or "residual AHC under zero mag- 
netization" reported in Pr2lr207 below w IK, if we 
assume small magnetic order of Ir 5d-electrons at . In 
fact, the AHC under the Ir spin structure in Fig. [TBI fc). 
which would be realized below under weak exchange 
field from "2in 2out" Pr order, is finite as shown in Fig. 
[TBI (a). The total AHC will be insensitive to the forma- 
tion of domain structure with three "2in 2out" Pr orders 
in Fig. [151 since (Tah's due to three [2 4, 2 f] structures 
are equivalent. The AHC obtained in the present study 
is expected to give a major part of the AHC observed in 
three dimensional compounds, as discussed in Appendix 
C. 

Since cr^'^ in the present model is nonzero unless 
'T'yi||'T's||''T'Cj the realization condition for the orbital 
mechanism is just the "non-coUinearity of the spin struc- 
ture", which is much more general than that for cr'^Aii- 
The orbital mechanism might be the origin of interest- 
ing spin structure-driven AHE in Fe3Sn2 [H, and 
PdCr02 
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Appendix A: Hopping integral between the sites 
with the different coordinates 



In this Appendix, we derive the hopping integrals be- 
tween the sites with the different d-orbital coordinates as 
shown in Fig. [S] Here, we represent the five d-orbitals 
xy, yz, zx, — and — as 1, 2, 3, 4 and 5. The 
wavefunctions in the c?-orbitals are given by 



V2^^ ^ ' 



-1 



1 

71 



2^ H 



<P5 = 

L r~ 

where 1^™ is the spherical harmonics and A — ^\hlA.-K. 

We consider the coordinate transformation matrix 
O'^^, which transforms {nx,ny, 71^)3 in the {xyz)^- 
coordinate into {nx,ny,nz)A in the (x?/2:)A-coordinate as 
{nx,ny,nz)Ad^^ = (n^;, ri^)B. It is given by 



-4 



-V 



qAB = 1 I 7 _4 4 



9 



(Al) 



-4 -1 



Since O^*-" and are equivalent to O^^, we have to 
derive only the hopping integral between sites A and B. 



Using rji 



^O^P where /, = x, y, z, the wavefunc- 



tion for orbital P at site B can be expressed as linear 
combination of the wavefunction for orbital 7 at site A. 
Thus, 



= ^aAB(/3,7)'/' 



where 



aAB(l, 1) 
aAB(l,2) 
aAB(l, 3) 
aAB(l,4) 

aAB(l, 5) 



— ^11 ^23 ' 



A7) 



^AB^AB 

^AB^AB 
"-^13 ^21 ) 

Of3^0,^2^ 



(A2) 



AB^AB 
13 ^^23 ' 



aAB(2, 1 
aAB(2,2 

aAB(2, 3 
aAB(2,4 

aAB(2, 5 

aAB(3, 1 
aAB(3,2 
aAB(3, 3 
aAB(3,4 

aAB(3, 5 



/^AB/--)Aa I /^AB/--)AB 
"^21 ^32 1 ^^22 ^31 
(^AB(^AB 
^^22 "-^33 

0^1^03^3" 

20^1^0^!^ + 

V^02^3^03^3^ 



(^AB(^AB 
^^23 "-^32 I 

0^3^03^^ 



= Ot^O^r 4 

= ^^03^3" 4 

_ r)AB^AB , 

— ^11 ^33 ^ 

= 730^3^0^3^. 



r)AB(oAB 
(^AB(^AB 
^AB^AB 

1- Ly^s 1^33 , 



Therefore, the hopping integral tsp.AaiRAB) = 
(B/3, RB\Ho\Aa, Ra) is given by 

tBl3,Aa{RAB) = ^ OAB 7)^A7,Aq (-Rab), (A3) 



where tA7,AQ(-RAB) = (A7, i^Bl^^olAa, i^A) is the usual 
hopping integral between the equivalent coordinates, 
which is given by the SK parameter table in Ref. (2^ . 



Appendix B: Local effective field from Pr 
tetrahedron 



In this Appendix, we derive the local effective field at 
Ir sites induced by the spin structure of Pr tetrahedron. 
Sites A, B, C and D of Pr tetrahedron are located at 
(1/4,0,0), (0,1/4,0), (0,0,1/4) and (1/4,1/4,1/4), re- 
spectively, in the xyz-coordinate as shown in Table I, and 
the center of the tetrahedron is located at (1/8, 1/8, 1/8). 

Under the strong field along [111](» O.TTesla), the 
spins of Pr tetrahedron form "3in lout" structure. The 
spin configurations at A^D are given by 



>34.iT 



„3m 



p3ilt 



(-l,l,l)/%/3, 
(1,-1,1)/V3, 
(1,1,-1)/V3, 
(1,1,1)/V3. 



In the intermediate field(^ O.TTesla), the spins of Pr 
tetrahedron is expected to form three kinds of "2in 2out" 
structures which have negative Zeeman energy. First, we 
consider the case in which only one of three "2in 2out" 
structures is realized. We choose one of three "2in 2out" 
structures, which is obtained by inverting only in 
the "3in lout" structure. That is, the configuration of 
the Pr spins in this "2in 2out" structure is given by 



,24.2t 
' A 
p242t 



,3ilt 
A ; 
,3ilt 



,24,2t 



,34,lt 



p2i2t 



-P^^'^ and 



Pq . We also consider another case where 
three "2in 2out" structure are averaged. Then, the Pr 
moments at sites A, B and C are given by 1/3 of the 
Pr moments in "3in lout" structure. That is, the Pr 
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spin in the averaged "2in 2out" structures is given by 



-'a /'^' -'^i 



B 



C 



73 



and P^^'^ = P^^'\ 

The effective magnetic fields at Ir sites are obtained by 
summing six Pr spins: Ir^ atom is surrounded by two Pr^, 
two Prfe and two Pr;, where we represent {i, j, k, 1} as a 
permutation of sites {A, B,C,D}. Therefore, the local 
exchange fields at Ir sites are given by 



(Bl) 



where Jdf is d-f exchange interaction. Hereafter, we as- 
sume Jdf > 0. We calculate the local exchange field h 
using above equation. The obtained results in each Pr 
structure are as follows: In the "Sin lout" structure case. 



h 



A = -J(3,l,l), 
hl^'^ = -J(l,3,l), 
-J(l,l,3). 



h 



34.1T 



^J(l,l,l). 



where J = 2Jcif/V3- In the "2in 2out" structure case. 



"■A 



h 



24.2t 
C 

24,2t 



-J(l,-1,3), 
-J(-l,l,3), 
-J(l,l,3). 
-J(-l,-l,3). 



In the averaged "2in 2out" structure case. 



"-A 



24.2t 
C 

2j2t 



-J(5/3,l,l), 
-J(l, 5/3,1), 
-J(l, 1,5/3). 
-J(l,l,l)/3. 



Next, we rewrite the obtained h's in the XyZ-coordinate 
using [nx^ny^nz] = {nx,ny,nz)0~^ , where the trans- 
formation matrix O^^ is given by 



o- 



1 



^-V3 -1 



-^1 -1 ^ 



(B2) 



In the "3in lout" structure case. 



h 



34.lt 



-J[-2V3,-2,5V2]/\/6, 
h^^^ = -J[2V3,-2,5\/2]/y6, 



34,lt 



-J[0,4,5V2]/V6, 



h: 



h^^^ = -J[0,0,3V2]/V6. 



In spherical coordinates, the direction of the local ex- 
change fields at site A, B and C are 9 = 29.5° and 



(0A,'/>B,</'c) = (-57r/6, -7r/6,7r/2). In the "2in 2out" 
structure case, 

h^j^'^^ = -J[-2\/3,6,3\/2]/\/6, 

hl^^^ = -J[2V3,6,3V2]/y6, 

h'^(j^'^ = -J[0,4,5V2]/V6, 

h^j^^^ = -J[0,8,V2]/y6, 

that is, 0a = 9b = 31.5°, Oc = 58.5°, and {cpA, '/'b, 0c) = 
(27r/3, 7r/3, 7r/2). In the average "2in 2out" structure 
case. 



I,2i2t 
"-A 

"-c 

34.1T 



h 



-J[-2V3, -2, ll\/2]/3\/6, 
-J[2\/3,-2,ll\/2]/3V6, 
-J[0,4, llV2]/3\/6, 
-J[0,0,%/2]/\/6, 



that is. 



14.4° and 



, '/>B,<^c) = 



(-57r/6,-7r/6,7r/2). 



Appendix C: AHC in three dimensional compounds 

In this paper, we have studied the AHE in the kagome 
lattice model, which represents the two-dimensional Mo 
or Ir network in the pyrochlore compounds. In the pres- 
ence of "3in lout" or "2in 2out" of Nd or Pr spin-ice 
order, it was shown that prominent spin structure-driven 
AHE are induced on the kagome lattice on the [1,1,1] 
plane. However, other three kagome layers on the [1,1,-1], 
[1,-1,1] and [-1,1,1] planes, which are not perpendicular 
to the magnetic field, also give finite contribution to the 
AHC. 

In this section, we shortly discuss the total AHC 
induced by four kagome lattices, assuming that these 
lattices are independent. Here, we put the magnetic 
field parallel to the [1,1,1] plane, which is given by 
the ABC plane in Fig. U] (a), and apply the electric 
field along Y axis. Then, the AHC due to the [1,1,1] 

plane, Cah'^'' given in the present study. Consid- 
ering the relative angles and positions of other three 
kagome lattices, it is easy to show that the total AHC 
due to the electric field on the [1,1,1] plane is given by 

<A - 4^H^'^' + + + ^H^'^'^Va. Note 

that [1, 1, —1], [1, —1, 1], and [—1, 1, 1] planes are respec- 
tively given by ABD, ACD, and BCD planes in Fig. H 
(a), where D represents the apical Ir site. 

First, we consider the AHC in Nd2Mo207. Under this 
120° structure of Ir spin (see Fig. [T5)) . the effective mag- 
netic flux due to the orbital AB effect is proportional to 
the tilting angle 6; ^I^'i^il a9. However, the orbital 
AB phase for other kagome layers are different. In the 
spinel- type kagome lattice studied in Ref. [2l|, we can 
show that the effective magnetic fiux for other layers are 
$[14,-1] = $[1,-1,1] = $[-1,1,1] ^ _(a/3)6l. Then, the 

total AHC is given by a"^^ = (2/3)cr]^H^'^'. 
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Next, we consider the AHC of Pr2lr207 under [2 I 
2 t]_ff-s-+o spin structure, considering only the con- 
ventional KL term that is proportional to the perpen- 
dicular magnetization. Then, it is easy to show that 

Thus, the 



^AH 



Jl,l,-1] 
'^AH 



-1,1,1] 



total AHC is given by a*i^^ = (4/3)ct 



'AH 

:i,i,i] 



AH 



As a result, a 



:i,i,i] 

AH 



gives the main contribution to the 
total AHC in both cases. Therefore, we expect that 
^ah'^' studied in the present study represents a major 
part of the AHC observed in three dimensional com- 
pounds. 
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